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Experiments by: H. Mabuchi and co-workers 
Norman Bridge Laboratory of Physics 12-33, California Institute of Technology, Pasadena, CA 91125, USA 

Quantum-limited estimation of an optical phase using adaptive (i.e. real-time feedback) techniques 
is reviewed. One case is explored in detail, as it can be understood using only elementary concepts 
such as photonic shot-noise and error analysis. Very recent experimental results are discussed. 

PACS numbers: 



I. INTRODUCTION 

Having been working on the topic of quantum-limit 
adaptive estimation of optical phase for some years now 
§ |, % 1 |, § 0, 1 1 , when I (HMW) was asked to con- 
tribute an article on it, the question "why now?" natu- 
rally arose. I believe there are three good answers to this 
question. In order of increasing importance, 

1. The work done by me, and more particularly my 
students, covers a wide range of cases. It is only 
now that it is possible to put them all in an overall 
context. 

2. Some of the latest results are actually the easiest 
to explain to a non-specialist. 

3. The first experiments verifying the theory have very 
recently been performed in the laboratory of Hideo 
Mabuchi at CalTcch. 

Before discussing these interesting developments, I 
should explain what quantum-limited adaptive phase es- 
timation means. The phase in question is an unknown 
(and possibly varying) optical phase ip. The aim is to 
estimate this phase, on the basis of measurements, as 
well as quantum mechanics allows. The allowed resources 
are practical ones: photodetectors, electronics, and linear 
(electro-)optics. The fundamental limitation is the num- 
ber of photons (per pulse or per coherence time) . There 
are also many practical limitations such as efficiency, the 
source of light, time delays, and processing limitations. 
These have all been considered || g, [j], but for sim- 
plicity I will not discuss them here. 

Where adaptive estimation comes in is that to realise 
the above aim using the allowed resources it turns out 
that it is necessary to use real-time feedback during the 
measurement. The basic idea is as follows. If ip is known 
to a very good approximation, then a simple measure- 
ment scheme will usually give near-optimal results. An 
example is homo dyne measurement of the phase quadra- 
ture of a pulse E3|. However if tp is totally unknown 



then the standard solution is to measure both quadra- 
tures (e.g. by heterodyne detection). This is far from 
optimal. The alternative is adaptive detection: use the 
results so far in the measurement to make an estimate ip 
of ip, and use this in a feedback loop to make the rest of 
the measurement closer to optimal. 

I believe this work is important for a number of rea- 
sons: 

First, quantum phase has a long and controversial his- 
tory 1 1 1 . Although ideal phase measurements can be 
defined, there is no way to make them without optical 
materials with arbitrarily high orders of nonlinearity Q . 
Hence it is of fundamental interest to know how well one 
can do with linear devices. 

Second, with continued miniaturisation of devices 
there will come a time when quantum limits do set the 
fundamental limits for technology. There is every reason 
to expect that phase estimation, as required for phase 
locking loops for example, will still be a practical con- 
cern then. 

Third, the idea of adaptive measurements may have 
applications far beyond that of phase measurements, so 
the theoretical and experimental expertise gained here 
potentially opens new frontiers in quantum measurement 
theory. It also has implications for the way in which 
we understand phenomena such as the "collapse of the 
wavef unction" . 

Finally, the cutting-edge technology required for 
quantum-limited adaptive measurements is a stepping 
stone towards more general methods for controlling quan- 
tum systems, another area of future importance as engi- 
neering heads towards the quantum realm. As will be 
discussed in Sec. V, the new technology being applied is 
both electronic and optical. 



II. OVERVIEW OF PAST RESULTS 

There are many different cases that can be considered, 
and they can be systematised by considering the follow- 
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ing four questions: 

1. Is the detection dyne (that is, using an optical local 
oscillator) or interferometric! 

2. Is the light source coherent (e.g. a laser), or non- 
classical (e.g. squeezed)? 

3. Is the scheme non-adaptive, or adaptive? 

4. Is it single-shot with a single phase, or CW with a 
varying phase? 

Two of these distinctions definitely require further expla- 
nation. 

First, dyne versus interferometric detection. These are 
illustrated in the figures below. 
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beams into the Mach-Zehnder interferometer are treated 
as quantum, and the photon number in both is counted 
H [§]. There is also a practical distinction in that dyne 
detection uses a photoreceiver that yields a photocurrent 
which does not distinguish individual photons, whereas 
interferometric detection uses photon counters. In both 
cases the aim is to measure the phase (p; in the former 
this phase is defined relative to the local oscillator, while 
for the latter it is as a relative phase between the two 
arms of the interferometer. 

Second, single-shot versus CW. In a single-shot mea- 
surement, there is a single pulse of light (which may 
however be split over two modes for the case of inter- 
ferometric detection). There is a single unknown phase 
ip imprinted on the pulse. The relevant parameter for the 
fundamental limit is n (or n), the (mean) number of pho- 
tons in the pulse. In Continuous- Wave detection, there 
is a beam (or two beams in the interferometric case) of 
light. A time-varying phase tp(t) is imprinted on (one of) 
the beam(s). In this article I will consider only the case 
where the time variation is that of white noise. That is, 



(p = y/K£(t). 



(2.1) 



This could arise from thermal mechanical fluctuations of 
miniaturised optical elements, for example. The param- 
eter k is the rate of phase diffusion, or, equivalently, the 
resulting linewidth of the beam. The relevant parameter 
in this case is N, the mean number of photons per coher- 
ence time: N = P/Hlvk, where P is the beam power. 

Having explained the various cases, I can now give the 
promised overview. The simplest results arise by con- 
sidering the mean square error (MSE) in the estimated 
phase, in the asymptotic limit where n or n or N go to 
oo. The results for the MSE scale as in table 1. 
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FIG. 1: Adaptive phase estimation: dyne detection (top) and 
interferometric detection (bottom). BS = Beam Splitter, and 
dashed lines are optical paths. 



TABLE I: Asymptotic mean square errors for phase estima- 
tion. The results that beat the standard quantum limits are 
underlined. 



It might well be objected that dyne detection (of 
which homodyne and heterodyne are two examples Hj) 
is merely one sort of interferometric measurement [|l2| , in 
that the local oscillator is interfered with the signal light 
before detection. This is quite true; I am using the two 
terms to make a distinction between how the accounting 
is done. In dyne detection only the signal pulse or beam 
is treated as quantum, and only the number of photons in 
it is counted for the purposes of determining the quantum 
limit. In interferometric detection both input pulses or 



There is a good deal of regularity in this table, as the 
reader may discern. I wish to draw attention to one fea- 
ture in particular. The sixteen MSEs (some still unde- 
termined) can be divided into four squares of four MSEs. 
In each square for which results are fully known, three 
of the MSEs scale in the same way. This scaling repre- 
sents the standard quantum limit (SQL) for that particu- 
lar case. The fourth, underlined in the above table, beats 
the SQL. In each case, this requires both non-classical 
light and adaptive detection. 
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III. LATEST RESULTS: CW DYNE DETECTION 

To try to explain adaptive phase estimation in more 
detail, I will concentrate now on the particular case of 
CW dyne detection. This is one the latest areas to be 
investigated, by Dominic Berry and myself but turns 
out to be probably the simplest to explain. To reiterate 
the basic idea, we want the best estimate of the current 
value of(p(t) which obeys ip = \^£,(t), where is white 
noise with unit spectral power. 

Consider first the case of coherent light of power P = 
frioa 2 , detected by interfering with a local oscillator at 
a balanced photoreceiver. The resulting dyne photocur- 
rent, suitably scaled, is Q 



I(t) = 2acos[$0) - <p(*)] + ((t), 



(3.1) 



where $(f) is the local oscillator phase, and ((t) is an- 
other Gaussian white noise term, independent of £(t), 
with unit spectral power. It can be thought of as local 
oscillator shot noise, or vacuum noise. In any case, it is 
quantum noise. 

The standard non-adaptive technique to estimate 
phase is to vary <E>(i) over all phases. This can be achieved 
by heterodyne detection with a detuning A 3> y/ctK, 
which makes $ vary as $(f) = <I>(0) + A x t. From 
Eq. (3.1), we would expect better sensitivity if we were 
to choose $(£) to maximise the slope of cos[$(i) — <p(t)]. 
That is, we should set $(t) = (p(t) + it/2 so that 



I(t) = 2asm[<f(t)-<f{t)]+((t) 
~ 2a[<p{t)-tp{t)]+£(t), 



(3.2) 



for (f(t) ~ <f{t)- The question is, how should we choose 
the estimate (p{t)l 

One obvious possibility is to choose it from the immedi- 
ately preceding photocurrent. We can rearrange Eq. (3.2) 
to get 



¥>(*) 



2a 



(3.3) 



For a given <p(t), we could thus form 
r t+6t 



fiimm(t + St) = — 



1 

St j, 



ds. 



(3.4) 



This has the MSE 



Vimmit + St) 



([ftmm(* + St) - tp(t + St)} 2 ) 
1 



Aa 2 Sf 



(3.5) 



This diverges as St — > dt, so clearly (^i mm is not a good 
estimate. Instead, we need a if>{£) that involves a finite 
time average. 

The optimal time-average for ip(t) can be determined 
as follows. Say the MSE in (p(t) is cr 2 {t). Over the interval 



[t, t + St), the diffusion of <p(t) causes this to increase to 

(*+**) = (m)-<p(t+st)] 2 ) 

= (mt)-v(t)] 2 ) + (Mt)-rtt+st)} 2 ) 

= a 2 (t) + KSt. (3.6) 

By standard error analysis, the optimal <p(t + St) weights 
ip(t) and (f>i mm (t + St) appropriately: 



ip(t + St) = a 2 (t + St) 
where 

1 



0imm(i + St) 



<p{t) 



(t + St) o- 2 old (t + St) 



a 2 (t + St) 



At + St) o- 2 old (t + St) 



(3.7) 
(3.8) 



Taking St — * dt yields a differential equation for a 2 (t) 
that has the stationary solution 



a 2 = 1/2VN , where N = a 2 / n. 



(3.9) 



Substituting this into Eq. (3.7) gives d(p{t) = (n/v 2 ) x 
(I(t)dt/2a). Since the feedback sets = <p(t) + tt/2, 
the feedback algorithm is simply 



$ = 



K I(t) 



(3.10) 



By way of comparison, for heterodyne (nonadaptive) 
detection, the stationary MSE can be shown to be 



a 2 et = 1/V2N. 



(3.11) 



That is, the adaptive technique offers a factor of 1/V2 
improvement in the mean square error. 



IV. WHAT IS QUANTUM ABOUT IT? 

This is a question I am often asked. A feedback algo- 
rithm of the form 



* = xHt)/2a 



(4.1) 



looks like a standard low frequency classical phase- 
locking algorithm with gain \. The quantum feature is 
that the optimal value of \ is finite: 



Xopt — 9 



2^/na. 



(4.2) 



In the classical limit (with no other noise sources), Xopt 
would be infinite, as can be seen by writing it explicitly 
in terms of classical quantities and h: 



Xopt = 2^/KP/hw. (4.3) 
For general x> the stationary MSE can be shown to be 



8a 2 2x 



(4.4) 
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which diverges as \ ~~ * 00 • Moreover, if \ differed from 
Xopt by a factor greater than about 2.4 (in either direc- 
tion), then all advantage gained by doing an adaptive 
detection would be lost. That is, the MSE would become 
greater than that which can be achieved by th e nonadap- 
tive technique of heterodyne detection ( 3.11 ). 

The advantage offered by adaptive detection is only 
a constant factor (of 1/V2) because so far I have dis- 
cussed only coherent light. To break the SQL scaling 
a 2 ~ N^ 1 / 2 it is necessary to use nonclassical light. The 
most obvious sort of nonclassical light to consider for 
this situation is broad-band squeezed light jOj. This is 
light where the spectrum of the dyne photocurrent of 
one quadrature, normalised as in Eq. ([Tl]), has a noise 
spectrum S below unity over some bandwidth broad com- 
pared to yftux. For moderate phase quadrature squeezing 
of spectral noise S, we can use the same feedback algo- 
rithm (4.1) to get an improvement in the MSE by a factor 
of y/S: 



a 2 = sfSjlsTN , for x = k/ct 2 



(4.5) 



It would appear from this that the MSE in the phase 
estimate could be made as small as one desires by making 
the phase quadrature more squeezed. However this is not 
the case because the noise gets "squeezed" into the am- 
plitude quadrature of the light jlO| . Because the lo- 
cal oscillator phase $ is only orthogonal to the estimated 
system phase tp, not the true phase ip, the amplitude 
fluctuations do contribute to the noise in the phase mea- 
surement. It can be shown that as a consequence there 
is an optimal non-zero squeezing spectrum S ~ iV -1 / 3 . 
This is, in a sense, a "more quantum" feature than the 
existence of shot noise which gave the above limit (3.9), 



because it is a consequence of the uncertainty relation 
between amplitude and phase for a light beam fl(i|) . In 
any case, with squeezed light we can beat the SQL, and 
find a new quantum limit for CW phase locking: 

a 2 ~ N- 2/z . (4.6) 



V. EXPERIMENTAL PROGRESS 

Performing a quantum-limited adaptive phase mea- 
surement has been one of the primary research objectives 
of Assoc. Prof. Hideo Mabuchi since his appointment at 
CalTech in 1999. As mentioned in the introduction, this 
is the first step towards a longer-term goal of investigat- 
ing quantum-limited control in general, and in particular 
in cavity QED systems. This first step has proven quite 
challenging in its own right! 

Rather than a CW measurement of a varying phase as 
discussed above, Mabuchi's group has been working on 
single-shot estimates of an unknown phase imprinted on 
a pulse of light, as originally considered in Ref. Q. Here 
the parameter n (see table 1) is simply the mean number 
of photons in a single pulse on which a single unknown 
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FIG. 2: Numerical results for the mean square phase error 
using optimal squeezing, for heterodyne (+) and adaptive (x) 
detection. Note that these were obtained using a more sophis- 
ticated feedback algorithm than that explained above. The 
analytical asymptotic results are: heterodyne ( — ), 0.66/iV 1 '' 2 
and adaptive (— ), 0.63/7V 2/3 . 



phase is imprinted. For an adaptive measurement, the 
phase must be estimated when only a fraction of these n 
photons have entered the detector. The fact that this can 
work even for n of order unity is described by Mabuchi 
as a "pretty wild fact that makes you think twice about 
field quantisation" . 

In this single-shot case the feedback algorithm is con- 
siderably more complicated than that in Sec. IV, as the 
local oscillator phase $(t) may depend nonlinear ly upon 
various integrals of the photocurrent during the course of 
a single pulse J(|. The form of these integrals can only 
be derived from the full quantum theory ||. Mabuchi 
thus decided to use digital electronics, and to discretise 
the feedback algorithm. 

Just as the CW result requires the magnitude of the 
feedback x to be close to its theoretical optimum (as dis- 
cussed in Sec. IV), the single shot case is sensitive to the 
design of the feedback loop. Through numerical model- 
ing, Mabuchi's group have found that the discretisation 
requires at least of order 100 divisions of the total pulse 
length T. 

The only electronics fast enough to do the required pro- 
cessing are field programmable gate arrays These 
are so fast (10s of MHz) that the feedback bandwidth, 
of around 5 MHz, is actually limited by the the band- 
width of the electro-optics (in particular the synthesiser 
used to modulate the feedback). This bandwidth allows 
each time-division to be of order 500ns, for a total time 
T of only about 50/xs. This still required the production 
of a coherent state of this duration. In other words, it 
required light with a shot-noise limited amplitude spec- 
trum down to about 50 kHz. This was another enormous 
technical challenge, as this is far below the usual fre- 
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quency regime of quantum-limited experiments UM . It 
was achieved using a very high Q cavity of lifetime 16 /is 
as a mode cleaner. 

Because they are not (as yet) using any non-classical 
source, the experimental signature Mabuchi's group is 
seeking is a reduction of the MSE (by at most a factor 
of 2 — see table 1) below the heterodyne limit of 0.5 /n. 
In work just released to the physics archive Pq| , this re- 
duction has been convincingly demonstrated. Of course, 
much work remains to be done. For example, Mabuchi 
plans to investigate the tails of the distribution as well 
as its standard error, to verify some of the predictions of 
Ref. @. Employing non-classical light to break the SQL 
is another longer term goal. Other experimental groups 
are interested in this as well. 



VI. CONCLUSIONS 

Estimating an unknown phase ip at the quantum limit 
is a difficult problem. There are many different cases that 



can be considered, but for all of them the standard quan- 
tum limits (SQLs) for the mean-square error arise for co- 
herent light or non-adaptive measurements. The SQLs 
can be beaten only by using nonclassical light in con- 
junction with adaptive measurements. When devices in- 
cluding phase-locking loops are sufficiently miniaturised, 
quantum limits will become practical limits. Thus, as 
well as being of theoretical interest, adaptive measure- 
ments should become part of the tool-kit of the future 
"quantum engineer" who seeks to manipulate quantum 
systems as well as nature allows. 

The improvement offered by adaptive detection over 
non-adaptive detection in phase estimation has very re- 
cently been achieved by the group of Hideo Mabuchi at 
CalTech. This experiment breaks new ground in the use 
of low-frequency optical mode cleaning and high-speed 
digital electronics in quantum- limited experiments. The 
theoretical and experimental techniques developed also 
pave the way for future research into the feedback-control 
of quantum optical systems in general. 
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